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Abstract 

We introduce a hierarchy of problems between the Dominating Set problem and the Power 
Dominating Set (PDS) problem called the f-round power dominating set (Ground PDS, for 
short) problem. For I = 1, this is the Dominating Set problem, and for I > n — 1, this is the 
PDS problem; here n denotes the number of nodes in the input graph. In PDS the goal is to 
find a minimum size set of nodes S that power dominates all the nodes, where a node v is power 
dominated if (1) v is in S or it has a neighbor in S, or (2) v has a neighbor u such that u and all of 
its neighbors except v are power dominated. Note that rule (1) is the same as for the Dominating 
Set problem, and that rule (2) is a type of propagation rule that applies iteratively. The ^-round 
PDS problem has the same set of rules as PDS, except we apply rule (2) in "parallel" in at most 
£ — 1 rounds. We prove that Ground PDS cannot be approximated better than 2 log ™ even for 
£ = 4 in general graphs. We provide a dynamic programming algorithm to solve £-round PDS 
optimally in polynomial time on graphs of bounded tree-width. We present a PTAS (polynomial 
time approximation scheme) for £-round PDS on planar graphs for £ = 0( ^°f^ n ). Finally, we give 
integer programming formulations for £-round PDS. 



1 Introduction 



The Power Dominating Set problem (PDS, for short) is a covering problem in which the goal is 
to "power dominate" (cover) all the nodes of a given undirected graph G by picking as few nodes as 
possible. There are two rules for power dominating the nodes; the first one has a "local" effect but 
the second one allows a type of indirect propagation. More precisely, given a set of nodes S, the set 
of nodes that are power dominated by S, denoted V(S), is obtained as follows 

(Rl) if node v is in S, then v and all of its neighbors are in V(S); 

(R2) (propagation) if node v is in V(S), one of its neighbors w is not in V(S), and all other neighbors 
of v are in V(S), then w is inserted into V(S). 

The PDS problem is to find a node-set S of minimum size that power dominates all the nodes (i.e., find 
Scy with \S\ minimum such that V(S) = V). For example, consider the planar graph G = (V,E) 
in Figure 1. The graph is obtained by taking the union of m paths on k + 1 nodes that all meet at a 
common node v; note that |V| = k ■ m + 1. It is easy to check that the size of a minimum PDS in G is 
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1; by taking S = {v} we get V{S) = V . In more detail, by applying (Rl) we power dominate v and 
all of its neighbors in the set X\\ after that, each node in the set X\ has exactly one neighbor that 
is not power dominated yet, namely, its neighbor in X^; thus, we can sequentially apply propagation 
rule (R2) to each node in X\ (in any order) to power dominate all the nodes in X<i; continuing in this 
way, all the nodes will be power dominated eventually. 

The PDS problem arose in the context of monitoring electric power networks. A power network 
contains a set of nodes and a set of edges connecting the nodes. A power network also contains a set of 
generators, which supply power, and a set of loads, where the power is directed to. In order to monitor 
a power network we need to measure all the state variables of the network by placing measurement 
devices. A Phasor Measurement Unit (PMU) is a measurement device placed on a node that has the 
ability to measure the voltage of the node and the current phase of the edges connected to the node; 
PMUs are expensive devices. The goal is to install the minimum number of PMUs such that the 
whole system can be monitored. These units have the capability of monitoring remote elements via 
propagation (as in Rule 2); see Brueni [1], Baldwin et al. [2], and Mili et al. [3]. Most measurement 
systems require one measurement device per node, but this does not apply to PMUs; hence, PMUs 
give a huge advantage. To see this in more detail consider a power network G = (V,E), and assume 
that the resistances of the edges in the power network are known, and the goal is to measure the 
voltages of all nodes. For simplicity, assume that there are no generators and loads. By placing a 
PMU at node v we can measure the voltage of v and the electrical current on each edge incident to v. 
Next, by using Ohm's law we can compute the voltage of any node in the neighborhood of v (Rule 1). 
Now assume that the voltage on v and all of its neighbors except w is known. By applying ohm's law 
we can compute the current on the edges incident to v except {v,w}. Next by using Kirchhoff's law 
we compute the current on the edge {v,w}. Finally, applying Ohm's law on the edge {v,w} gives us 
the voltage of w (Rule 2). 

PMUs are used to monitor large system stability and to give warnings of system-wide failures. 
PMUs have become increasingly popular for monitoring power networks, and have been installed by 
several electricity companies since 1988 [4, 5]. For example, the USA Western States Coordinating 
Council (WSCC) had installed around 34 PMUs by 1999 [5]. By now, several hundred PMUs have 
been installed world wide [6], Some researchers in electrical engineering regard PMUs as the most 
important device for the future of power systems [7]. 

The PDS problem is NP-hard even when the input graph is bipartite [8]. For further references for 
the PDS problem please see [8, 9, 10, 11, 12, 13, 14, 15]. PDS is a generalization of the Dominating 
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Figure 1: 

Set problem. A Dominating Set of a graph G = (V, E) is a set of nodes S such that every node 
in the graph is either in S or has a neighbor in S. The problem of finding a Dominating Set of 
minimum size in a given graph G has been studied extensively in the past 20 years, see the books 
by Haynes et al. [16, 17]. The Dominating Set problem is a well-known NP-hard problem [18]. A 
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simple greedy algorithm achieves a logarithmic approximation guarantee, [19], and modulo the P 7^ NP 
conjecture, no polynomial time algorithm gives a better approximation guarantee, [20, 21, 22]. 

In this paper we introduce a hierarchy of problems between Dominating Set and PDS, by adding 
a parameter £ to PDS which restricts the number of "parallel" rounds of propagation that can be 
applied. The rules are the same as PDS, except we try to apply the propagation rule in parallel as 
much as possible. In the first round we apply the rule (Rl) to all the nodes in S, and for the rest of the 
rounds we only consider "parallel" application of the propagation rule (R2). In every "parallel" round 
we power dominate all the new nodes that can be power dominated by applying the propagation rule 
to all of the nodes that are power dominated in the previous "parallel" rounds. Given a parameter 
£, the £-round PDS problem is the problem in which we want to power dominate all of the nodes in 
at most £ parallel rounds. Clearly, the Ground PDS problem for £ = 1 is exactly the Dominating 
Set problem, and for a graph G with n nodes the Ground PDS problem for £ > n — 1 is exactly the 
PDS problem. The notion of parallel propagation comes from the fact that changes in the electrical 
network propagate in parallel and not sequentially. A feasible solution for the PDS problem provides 
a plan for installing monitoring devices to monitor the whole power network, but it does not provide 
any guarantees on the time-lag between a fault in the network and its detection. Deducing information 
through a parallel round of propagation takes one unit of time and in some applications we want to 
detect a failure in the network after at most I units of time. The addition of the parameter £ achieves 
this time constraint. 

The practical motivation for the Ground PDS problem has been explained above. In addition, there 
are some theoretical motivations. Although the PDS problem has been studied since 1993, there 
are very few algorithmic results (including approximation algorithms). The Ground PDS problem 
serves as a unified model for studying the PDS problem and the Dominating Set problem. The 
introduction of the parameter I allows us to examine the complexity of the problem in terms of £: how 
does the threshold for the hardness of approximation vary in terms of I ? The hardness threshold is 
logarithmic for £ = 1 (the Dominating Set problem) and it is f2(2 log '") for £ = n — 1 (the PDS 
problem). Is the latter hardness threshold valid for constant £1 Moreover, for planar graphs there is 
extensive recent literature on PTASs (polynomial time approximation schemes) for the Dominating 
Set problem and its variants, but these results do not apply to the PDS problem. A major open 
question in the area is whether there exists a PTAS for the PDS problem on planar graphs. One 
avenue that may lead to advances on this question is to design a PTAS for the Ground PDS problem 
on planar graphs, for small values of £. Integer programming formulations for the PDS problem have 
been studied, and we give a new formulation in the last part of this paper. Our formulation is based 
on the notion of parallel propagation in £ rounds. We first give a formulation for the f-round PDS 
problem, and then modify it to get another formulation for the PDS problem. 

1.1 Our main results 

We initiate the study of a natural extension of the PDS problem and prove the following main results. 

• For general graphs, we present a reduction from the MinRep problem to the ^-round PDS 
problem which shows that f-round PDS for £ > 4 cannot be approximated better than 2 log tn , 
unless NP C DJ\ME(n polyl09 ^). We use a reduction similar to one that has been used to prove 
the same hardness result for the PDS and the directed PDS problems in a paper jointly authored 
with M. Stilp [13]. 

• We provide a dynamic programming algorithm to solve the £-round PDS problem optimally in 
polynomial time on graphs of bounded tree-width. This dynamic programming algorithm is 
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based on our reformulation for the £-round PDS problem. This reformulation is an extension 
of the one introduced by Guo et al. [10] for PDS. Guo et al. [10] gave a new formulation for 
PDS in terms of "valid orientation" of the edges; they use this formulation to design a dynamic 
programming algorithm to solve the PDS problem optimally in linear time on graphs of bounded 
tree- width. 

• We focus on planar graphs, and give a PTAS for f-round PDS for £ = 0( lo 1 °^ n ). Baker's PTAS 
[23] for the Dominating Set problem on planar graphs is a special case of our result with £ = 1, 
and no similar result of this type was previously known for £ > 1 . We also show that the Ground 
PDS problem in planar graphs is NP-hard for all £ > 1. Note that our PTAS does not apply to 
PDS in general, because the running time is super-polynomial for £ = ^ dJ"^ ). 

• Finally we study integer programming formulations for ^-round PDS. 

Here is a brief discussion on the relation between this paper and some previous joint work with 
M. Stilp on the PDS problem [13]. Although the results in [13] do not imply any of the results in 
this paper, there are two topics that uses similar methods: (1) The hardness result for ^-round PDS 
(Theorem 5.1) extends the construction used to prove the hardness of directed PDS in [13], (2) The 
dynamic programming algorithm for ^-round PDS (Section 3) and the algorithm for directed PDS in 
[13] are both based on reformulation of the problems that extend the methods of Guo et al. [10]. 



2 Preliminaries 

Most of the graphs that we consider here are undirected. Given a graph G = (V,E), we denote the 
number of nodes in the graph by n. Sometimes the graphs that we consider have some directed edges 
in addition to undirected edges. Given such a graph, we denote by d~(v) and d + {v) the number of 
directed edges with v as the head and the tail respectively. Let G denote the underlying undirected 
graph obtained from G by ignoring the direction of the directed edges. Then the closed neighborhood 

of a node u in G is defined by N[u] = |i> : {u,v} € £7(Cr)| U {u}. Now we define the "parallel" 
propagation rule formally. 

Definition 2.1 Given a graph G = (V,E) and a subset of nodes SCF, the set of nodes that can be 
power dominated by applying at most k rounds of parallel propagation, denoted by V k (S), is defined 
recursively as follows: 

I W \ V k ~ l (S) |J [v : (u, v) 6 E, N[u] \ {v} C V k ~ 1 (S)} k>2 
We can now define the Ground PDS problem formally. 

Definition 2.2 (£-round PDS) Given a parameter £, the ^-round PDS problem is the problem in which 
we are given a graph G = (V, E) and the goal is to find a minimum size subset of nodes S C V , such 
that V e (S) = V. 

Given a graph G and a parameter £, we denote by Opt^(G) the size of the optimal solution for the 
Ground PDS problem. It is easy to see that the size of an optimal solution for ^-round PDS does 
not increase by increasing the value of the parameter I; to see this, consider an optimal solution S* 
for ^-round PDS and note that it is also a feasible solution for (£ + l)-round PDS. This proves the 
following property of the Ground PDS problem. 
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Proposition 2.3 Let G = (V, E) be a graph with n nodes, then for any parameter 1 <£<£'< n we 
have Opt e ,(G) < Opt e (G). 

The £-round PDS problem is very different from the Dominating Set problem, even for £ = 0(1). 
The following result shows that the size of an optimal solution for the Dominating Set problem can 
be much bigger than the size of an optimal solution for the 2-round PDS problem. 

Proposition 2.4 Let £ > 1 be a given parameter, the ratio between the size of the optimal solution 
for l-round PDS and {£ + l)-round PDS can be Q( j) even on planar graphs. 

Proof: Consider the graph G that is obtained by taking m paths, P±, . . . ,-P m , each of length £ + 1 
that share a common node v as shown in Figure 1 (with k = £ + 1). Note that G has n = m • (^+ 1) + 1 
nodes. It is easy to check that v can power dominate the entire graph in exactly £ + 1 parallel rounds; 
in the first round v and X\ are power dominated, in the second round Xi is power dominated, and so 
on. This shows that Opt^ +1 (G) = 1. Next we prove that the size of the optimal solution for ^-round 
PDS is at least m = jtj = O(j). For the sake of contradiction assume that it is not, so there is an 
optimal solution S* such that \S*\ < m. This implies that there is a path Pi that has no node in 
S* except possibly the node v. This means that the last node Ui on the path Pi (in Xk) cannot be 
power dominated since there are only £ parallel rounds and the distance between v and Ui is £ + 1. 
This contradicts the fact that S* is a feasible solution for ^-round PDS, so Opt^(G) > m. Also note 
that by taking S to be the set of middle nodes of the paths Pi, ... , P m we can power dominate all of 
the nodes in £ rounds (in fact, in | rounds). This shows that Opt^(G) = m, and therefore the ratio 
between Opt £ (G) and Opt e+1 (G) is m = 0(j). □ 

The f-round PDS problem can be generalized as follows. Given a subset V of the nodes, find a 
minimum size set of nodes that power dominates the nodes in V' in at most £ parallel rounds. We will 
use this generalization in our PTAS in Section 4. 

Definition 2.5 (Generalized ^-round PDS) An instance of the generalized ^-round PDS problem is 
given by a pair (G,V), where G = (V,E) is an undirected graph and V C V, and the goal is to find 
a minimum size set of nodes S such that V C V (S). 

3 Dynamic programming for f-round PDS 

In this section we provide a dynamic programming algorithm to solve the generalized Ground PDS 
problem optimally. Our dynamic programming algorithm is based on tree decompositions [24]. The 
running time of this algorithm is exponential both in the tree-width of the given graph and the 
logarithm of the parameter £. Therefore, we obtain a polynomial time algorithm for graphs of bounded 
tree-width. Our dynamic programming algorithm is based on our reformulation for the £-round PDS 
problem. This reformulation is an extension of a new formulation for PDS introduced by Guo et al. 
[10]. Our dynamic programming is similar to the dynamic programming for the Dominating Set 
problem given in [25] and the dynamic programming for PDS given in [10]. 

Definition 3.1 A tree decomposition of a graph G = (V,E) is a pair ({Xi C V \ i e /} , T), where 
T = (L, F) is a tree, satisfying the following properties: 

1. \J ia Xi = V; 
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2. For every edge {u, v} G E there exists an i G / such that {u,v} C Xi\ 

3. For all i,j, k G I, if j is on the unique path from i to k in T, then we have: Xi D X^ C Xj. 

TTie width 0/ ({Xj | i £ 1} ,T) is defined as max ie j\Xi\ — 1. TTie tree-width 0/ G is defined as the 
minimum width over all tree decompositions. The nodes ofT are called T- nodes, and each Xi is called 
a bag. 

For designing a dynamic programming algorithm based on tree decomposition it is usually easier to 
work with nice tree decompositions [26] that have a simple structure defined as follows. 

Definition 3.2 A tree decomposition {{Xi\i G I},T) is called a nice tree decomposition if the follow- 
ing conditions are satisfied: 

1. T is rooted, at r, and every node ofT has at most 2 children; 

2. If a node i G I has two children j and k, then Xi = Xj = Xj. (in this case i is called a Join 
Node,); 

3. If a node i has one child j, then either of the following holds: 

(a) Xj C Xi and |Xj \Xj\ = l (in this case i is called an Insert Node,). 

(b) Xi C Xj and \Xj \ Xj| = 1 (in this case i is called a Forget Node,). 

A tree decomposition can be transformed into a nice tree decomposition [26]. Therefore, to design 
a dynamic programming algorithm based on tree decomposition, we can assume that we are given a 
nice tree decomposition of the input graph. 

Lemma 3.3 ([26]) Given a tree decomposition of a graph G that has width k and O(n) nodes, where 
n is the number of nodes of G, a nice tree decomposition of G that also has width k and 0(n) nodes 
can be found in time 0(n). 

Now we provide a reformulation, called timed- orientation, for (generalized) £-round PDS in terms 
of the orientation of the edges and labeling of the nodes. This reformulation makes it possible to 
design a dynamic programming algorithm to solve the (generalized) Ground PDS problem optimally 
in polynomial time on graphs of bounded tree- width. In this formulation we orient the edges in order to 
show the direction that the propagation rule is applied; if the node w is power dominated by applying 
the propagation rule on v then we orient the edge {v, w} from v toward w. Moreover, if w is not power 
dominated through v we will leave the edge {v, w} undirected. In order to keep track of the round in 
which every node is power dominated, we introduce a time vector {t v : v G V}. The round in which 
v is power dominated is denoted by t v and it can take any value from the set {0, ...,£} U {+00}; the 
+00 is used to denote a node that will not be power dominated (this is needed for generalized Ground 
PDS). Consider a directed edge (it, v). There are two cases: either t v = 1 or t v > 1. In the former case 
u should be in the (optimal) solution, since t v = 1 means v is power dominated at the first round by 
applying the domination rule. In the latter case u and all of its neighbors except v (i.e N[u] \ v) should 
be power dominated before we can apply the propagation rule. In fact, v will be power dominated 
right after the round where the last node in N[u] \v is power dominated (see the property P5 in the 
following definition). 
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Definition 3.4 (valid timed-orientation) Let (G, V) be an instance of the generalized l-round PDS 
problem where G = (V, E) is a graph and V' C V is a subset of nodes. A valid timed-orientation for 
(G, V) is a graph G Q = (V, Ed U E u ) such that for every {u, v } G E either there is a directed edge 
(u,v) or (v,u) in Ed or an undirected edge {u, v} in E u , together with the time vector {t v : v G V} 
(possible values for t v are {0, 1, ■■■,£} U {+00} ) that satisfies the following properties: 



(PI) W G V : < t v < I, 

(P2) Vv G V : 1 < t v < £ => d~(v) = 1, 

(P3) Vv€V :t v = +oo=> dr(v) = d + (v) = 0, 

(P4) Vv G V : t v = cT(u) = 0, 



(P5) V(u, v) E Ed ■ t v 



1 */t„ = 

1 + maxjiy, : w G iV[u] — v} otherwise. 



The set O = {v G V : t v = 0} is called the origin of the valid timed-orientation. 



Now we show that the existence of a valid timed-orientation with S as the origin is equivalent to 
having S as a feasible solution to the ^-round PDS problem. 

Lemma 3.5 Let G = (V, E) be an undirected graph, and (G, V') be an instance of the generalized i- 
round PDS problem. Then S <^V is the origin of a valid timed- orientation if and only ifV C V £ (S). 

Proof: Suppose S C V power dominates V' in at most £ parallel rounds. Now we construct a valid 
timed-orientation with S as the origin. We orient the edges in the same way as the propagation rule 
applies. Consider an edge {u, v} and assume that v is power dominated by applying power domination 
rules to u. Then we orient the edge {u, v} from u toward v. Let Ed be the set of oriented edges and 
let E u be the rest of the edges. This defines the graph G Q = (V, Ed U E u ). Now we define the time 
vector. First, for any node v that is not power dominated, set t v = +00. Second, for every v G S, set 
t v = 0. Finally, define t v for the rest of the nodes as t v = min{r > 1 : v £ T >r {S)}. It is straightforward 
to check that the above orientation and the time vector satisfy all of the properties PI to P5 in the 
definition of the valid timed-orientation. Therefore, there is a valid timed-orientation with S as its 
origin. 

Now assume that G has a valid timed-orientation with S as the origin. So there is an oriented graph 
G = (V, EdL)E u ) and a time vector {t v : v G V} that satisfy the properties PI to P5 of the valid timed- 
orientation. Define V r = {v G V : < t v < r}. We now prove by induction that V r C V r {S) which 
implies that C V^(S). First note that this proves the lemma, since any node v G V has < t v < £ 
by the property PI. Note that any node v with t v = is in S, so it is immediately power dominated. 
Also any node v with t v = 1 is power dominated (by the domination rule) in the first round, since by 
P2 it has an incoming directed edge, say (w,v), and by P5 the node w is in S. Hence, V\ C V l {S). 
Now assume that the induction hypothesis holds for r = k (where k > 1), that is Vk C V k (S). We will 
show that it also holds for r = k + 1. Consider a node v such that t v = k + 1. By P2 it has exactly one 
incoming edge, say (u, v) G Ed- The node u has t u > 1 (since if t u = then P5 implies t v = 1), so by 
the property P5 any w G N[u] — v has t w < k. Therefore, N[u] \ {v} C V k (S). This is correct for any 
v with t v = k + 1, so X = {v G V : t v = k + 1} C {v : (u,v) G E,N[u] \ {v} C V k (S)}. Therefore, by 
the definition of the parallel propagation X can be power dominated in one parallel round, so we have 
V k+1 = V k UX CV k+1 {S). This proves the induction step and completes the proof. □ 
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Theorem 3.6 Given a pair (G, V) where G = (V, E) is a graph with tree-width k and V' C V , a 
minimum size set S C V such that V C V e (S) can be obtained in time O(c me+k ^ ogi ■ \V\), for some 
global constant c, where m e is the maximum number of edges that a bag can have. 

Proof: Please refer to Appendix A. 2 for the details of the dynamic programming. □ 

In general graphs the number of edges in a bag is at most ( k ~^ 1 ) but in planar graphs, it is linear 
in the number of nodes, where k denotes the tree-width. Therefore, the above theorem implies the 
following corollary. 

Corollary 3.7 Given a pair (G, V) where G = (V, E) is a planar graph with tree-width k and V C V , 
a minimum size set S C V such that V' C V (S) can be obtained in time 0(c fclog ^ • |V|), for some 
global constant c. 

4 f-round PDS on planar graphs 

In this section we present a PTAS (polynomial time approximation scheme 1 ) for the Ground PDS 
problem on planar graphs when I = 0( ^°f ( ^ n ). Baker's PTAS [23] for the Dominating Set problem 
on planar graphs is a special case of our result with £ = 1, but there are no previous results of this 
type for £ > 1. Our PTAS works in the same fashion as Baker's PTAS, but our analysis and proofs 
are novel contributions of this paper. Demaine and Hajiaghayi [27] recently used bidimensionality 
theory to design PTASs for some variants of the Dominating Set problem on planar graphs (such 
as Connected Dominating Set), but their methods do not apply to the Ground PDS problem 
because the relevant parameter is not bidimensional. We also have the following NP-hardness result 
for Ground PDS on planar graphs. 

Proposition 4.1 For any £ > 1 the l-round PDS problem is an NP-hard problem even on planar 
graphs. 

Proof: We use a modification of the reduction that has been used to prove that the PDS problem on 
planar graphs is NP-hard [9, 10]. Please refer to Appendix A.l for more details. □ 

Now we describe our PTAS for the £-round PDS problem on planar graphs, when £ is small. First 
we provide some useful definitions and notations. Consider an embedding of a planar graph G. We 
define the nodes at level i denoted by Lj as follows [23]. Let L\ be the set of nodes on the exterior 
face in the given embedding of G. For i > 1, the set Li is defined as the set of nodes on the exterior 
face of the graph induced on V \ Vf^Lj. We denote by L(a, b) = \J b i=a Li the set of nodes at levels a 
through b. A planar graph is called k-outerplanar if it has an embedding where no node is at level 
greater than k. For example, consider the graph in Figure 2. Clearly, the graph is a 2-outerplanar 
graph. The set L\ = {u±, U2, ■ ■ ■ , u&} is the set of nodes at level 1 and the set L2 = {i>i, V2, . . . ,vs} is 
the set of nodes at level 2. Given a graph G = (V, E) and V C V, we denote the subgraph induced 
on V by G[V']. 

Before describing our PTAS, let us look at Baker's PTAS [23] for the Dominating Set problem 
on planar graphs. Given the parameter e = jk Baker's algorithm finds a feasible solution with size 
within (1 + e) times the optimal value. This algorithm considers k different decompositions T*i, . . . ,T>k 

1 A polynomial time approximation scheme (PTAS) is an algorithm that given any fixed e > provides a solution with 
cost within + times the optimal value in polynomial time. 



s 



of the nodes of G and then finds a feasible solution for each of them. The ith decomposition consists 
of blocks with k + 1 consecutive levels. The jth. block in decomposition T>i contains nodes of levels 
jk + i through (j + l)k + i (note that each T>i is obtained from T>\ by shifting the levels). Also note 
that every two consecutive blocks in a given decomposition share a common level. Next, the algorithm 
solves the Dominating Set problem optimally on each block of £>j. This is possible since each block 
is (k + 1) -outer planar; hence, it has tree- width at most 3(/c + 1) — 1 and the Dominating Set problem 
can be solved optimally by dynamic programming [23] (the result in Section 3 for t = 1 also shows 
this fact). Then, it takes the union of the optimal solutions for the blocks in T>i to obtain a feasible 
solution for the Dominating Set problem in the graph G. Let Si denote this feasible solution. 
Finally, the algorithm outputs the solution that has minimum size, i.e. minj = i among all k 

decompositions. It is not hard to argue that the size of this solution is within = 1 + e times the 
optimal value. The key property that is needed for this argument is the fact that consecutive blocks 
share a common level. 



Now let us describe our PTAS informally. Consider a parameter k, that is a function of the pa- 
rameter I and the approximation factor (1 + e); k will be defined later in the formal description of 
our algorithm. We decompose the graph in k different ways, T>i, . . . , in each decomposition the 
graph is decomposed into blocks of k + 41 — 2 consecutive levels. The jth block in T>i is defined 
as B is j = L(jk + i - 21 + 1, (j + l)k + i - 1 + 21 - 1). We denote the k middle levels of B itj by 
C%a = L(Jk + i, (j + \)k + i — 1). In our PTAS, for each decomposition T>i, we optimally solve instances 
of the generalized Ground PDS problem (2jj = [G[Bij], Cij}) defined for each block in £>j. Note 
that each instance lij can be optimally solved by using the dynamic programming algorithm given 
in Section 3. Let Oij denote the optimal solution for this instance. Then we take the union of the 
solutions corresponding to blocks in T>i, LTj = UjxjCi.j- By doing this for all k decompositions we get 
k different feasible solutions for the original graph G. Finally, we choose the solution with minimum 
size among these k solutions. We will see that the 2^ — 1 extra levels around the k middle levels and 
the common levels between consecutive blocks plays an important role in the feasibility and the near 
optimality of the final output of the algorithm. 

To make the role of common levels clear, consider an instance of 4-round PDS shown in Figure 
2. Assume that we partition the planar graph into two levels. The first level is the outer cycle and 
the second level is the inner cycle. It is easy to check that the size of an optimal solution in any 
one of these partitions is 1. For example {u±} and {^5} are optimal solutions for partition 1 and 
partition 2 respectively. But if we consider the original graph G, it is straightforward to check that 
S = {ui} U {^5} is not a feasible solution for the instance G of 4-round PDS; the set of nodes that can 
be power dominated in at most four parallel rounds is V (S) = {tti, U2, U5, Us, v±, V4, V5, Vq}. Note that 
the propagation rule was applied in the subgraph but not in the original graph, so in the subgraph u 




Figure 2: A 2-outerplanar graph 
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may have all but one node of N[u] in V l (S) but this need not hold for the original graph (see us, for 
example) . 

To prevent such a problem, we need to consider extra levels around each block as we did in the 
blocks Bi j above. We will show that a feasible solution, found by our algorithm, for the instance Zj ~ 
will power dominate at least its k middle levels in the original graph G. This implies that the union of 
the solutions for the blocks in a given decomposition will be a feasible solution for G, since the union 
of the k middle levels of the blocks covers all the nodes in G. 

Algorithm 1 PTAS for Ground PDS 
1: Given a planar embedding of the graph G, and the parameter < e < 1. 
2: Let k = 4- f|]. 
3: for i = 1 to k do 
4: for all j > do 

5: Solve "generalized" Ground PDS on (G[B itj ], dj) 

6: Let Oij be an optimal solution for (G[Bij],Cij) 

7: end for 

8: 11; = ^j>oO id 

9: end for 

10: r <— ar grain {\Tli\ : i = 1, • • • , k} 

ii: Output n = n r . 



Theorem 4.2 Let £ be a given parameter, where £ = 0( lo °^ n ). Then Algorithm 1 is a PTAS for 
the £-round PDS problem on planar graphs. 

Proof: Let IT* be an optimal solution for Ground PDS in G. To prove the theorem, it is enough to 
prove the following two claims: 1) Ojj is a feasible solution for the instance {G, Cjj) of the generalized 
Ground PDS problem, 2) IT* n Bij is a feasible solution for (G[Bij],Cij). First let us see how the 
theorem follows from the above claims. The first claim shows that IL, for each i, is a feasible solution for 
the Ground PDS problem for G, since Uj>oCij = V. The second claim shows that \Oij\ < |H* n Bij\, 
so |ILJ < Y2j |n* H Bij\. In the right hand side we counted the nodes in the optimal solution twice 
on 41 — 2 common levels between any two consecutive blocks. By considering all values of parameter 
i, 1 < i < k, we can find an i such that the number of double counted nodes in the optimal solution 
IT is at most ^ |n*|. This implies that |H | < (1 + ^) |H* |, so by setting k = 4 • \f\ we get a 
(1 + e)-approximation algorithm. 

Now we analyze the running time of the algorithm. In step (5) we solve an instance of the generalized 
Ground PDS problem. The graph in this instance has k-\-A£ — 2 levels, so it is a {k + Al — 2)-outerplanar 
graph. It is known that the tree- width of any d-outerplanar graph is at most 3d — 1 [28]. Therefore step 
(5) of our algorithm can be done in c °(^/ elo s^) time by Corollary 3.7, since the graph has tree-width 
at most 3(k + 4£) < 12([|] +£). This shows that £ should be 0( Io 1 g 1 ^ n ) in order to have a polynomial 

time algorithm and in this case step (5) can be done in n°W time. Also note that the value of j can 
be at most j, since the number of levels of a planar graph is at most n (the number of nodes) and 
each Cij has k levels. This shows that the algorithm executes step (5) at most k x f = n times. Also 
the steps (8) and (10) can be done in polynomial time. Therefore the running time of the algorithm 
is polynomial time in the number of nodes of G for a fixed e. 

Now we prove the two claims stated above: 1) Oij is a feasible solution for (G, Ci,j), 2) IT* n Bij is 
a feasible solution for (G[Bij],Cij). 
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Proof of the first claim: We know that Oij is a feasible solution for (G[Bij], Cij). Let t v denote 
the round in which v was power dominated in G[Bij]. So any node v € Cij and possibly some nodes 
v G Bij \ Ci j satisfy < t v < I. For simplicity we use L s to denote the levels L(Jk + i — 2(£ — s), (j + 
l)k + i — 1 + 2(i — s)) for any s > 1, and also L° = Bij. Observe that L s+1 (for s > 1) is obtained 
from L s by deleting the first two levels of L s and the last two levels of L s . 

First note that by taking Oij all nodes v with t v = in L° are power dominated in the graph G. 
Now we claim that the following statement is correct: for each s, 1 < s < £, all nodes v with t v < s in 
L s are power dominated in the graph G. We prove the statement by induction on s. The base case 
s = 1 is trivial by applying the first rule of PDS, since any node v E L 1 with t v = 1 had a neighbor 
ii 6 L° with = 0. Therefore all nodes with t v < 1 in L 1 are power dominated. Assume that the 
statement is correct for all s < s'. Consider a node v with i„ = s' which lies in L s , and assume that 
it was power dominated by applying the propagation rule on u. It is easy to see that u is inside the 
levels L(jk + i — 2{i — s') — 1, (j + l)k + i — 1 + 2(£ — s') + 1), and therefore all neighbors of u are inside 
the levels L(jk + i - 2(£ - s') - 2, (j + l)fc + i - 1 + 2(£ - s') + 2) = L 5 '" 1 (see Figure 3). As v was 
power dominated by u, any zz; 6 N[u] — v satisfies t w < s' — 1. By the induction hypothesis any node 
w G JV[ti] — u is already power dominated in the graph G. The propagation rule can be applied to 
u to power dominate v. This completes the induction step and proves the statement. An important 
property that should be preserved is that all nodes v with t v = s' should be power dominated in 
parallel. It is easy to see that in the above proof any such node can be power dominated in parallel 
since they were power dominated in parallel in the induced graph G[Bi j]. 

To prove the first claim it is enough to note that the set L is exactly Cj.j, so all the nodes in Cij 
with < t v < £ (which is exactly all the nodes in Cij) can be power dominated in the graph G. 




I s '- 1 



Figure 3: Induction step 

Proof of the second claim: We know that II* is a feasible solution for G. Let t v denote the round 
in which node v G Bij was power dominated in G. So each node v G Bij satisfies: < t v < £. Define 
L s as before. The same induction hypothesis as above will prove the statement: for each s, < s < £, 
all nodes v with t v < s in L s can be power dominated in the induced subgraph G[Bi j]. The proof is 
similar to the first claim. Note that the set L is exactly Cij, so all the nodes in dj can be power 
dominated in G[Bij] in at most £ parallel rounds. □ 

5 Hardness of ^-round PDS 

In this section we prove the following result by a reduction from the MlNREP problem. The reduction 
given here is similar to the reduction used to prove the same hardness of approximation for the Directed 
PDS problem [13]. 
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Theorem 5.1 The i-round PDS problem for any £ > 4 cannot be approximated within 2 logl ratio, 
for any fixed e > 0, unless NP C DTIME{nP ol y lo 9^). 

In the MinRep problem we are given a bipartite graph G = (A,B,E) with a partitioning of 
j4 and B into equal size subsets, say A = U?=i ^« an d ^ = Ui=i^«) where \Ai\ = tua = j^j 
and \Bi\ = tub = This partitioning naturally defines a super bipartite graph TL = {A,B,£). 
The super nodes of TL are A = {A\, A2, • • • , A qA \ and B = {B\, B2, ■ ■ ■ , B qg }, and the super edges 
are £ = {AiBj\3a G Ai,b G -Bj : {a, b} € -E(G)}. We say that a super edge AjSj is covered by 
{a, b} G E(G) if a G and 6 G -Bj. The goal in MinRep is to pick the minimum number of nodes 
A' U B' C y(G) from G to cover all the super edges in TL. The following theorem states the hardness 
of the MinRep problem [29]. 

Theorem 5.2 [29] The MinRep problem cannot be approximated within the ratio 2 1091 tn , for any 
fixed e > 0, unless NP C DTIME(nP ol y lo 9^), where n = \V(G)\. 

The reduction: Theorem 5.1 is proved by a reduction from the MinRep problem. In the following 
we create an instance G = (V, E) of Ground PDS from a given instance G = (A, B, E)(TL = (A, £>, £)) 
of the MinRep problem. 

1. Add a new node w* (master node) to the graph G, and add an edge between w* and all the 
other nodes in G. Also add three new nodes w^w^Ws and connect them to it?*. 

2. Vi G {1, . . . , ?a} , j G {1, . . . , q B } do the following: 

(a) Let Eij = {e±, e?, ■ . ■ , e K } be the set of edges between Ai = {a^, . . . ,a,i mA } and Bj = 
{bj 1 , . . . , bj mB } in G, where k is the number of edges between Ai and Bj. 

(b) Remove E^ from G. 

(c) Let the edge e q G Eij be incident to ai q and bj q (in G). In this labeling for simplicity the 
same node might get different labels. Let Dij be the graph in Figure 4 (a dashed line shows 
an edge between a node and the master node w*). Make A = 4 new copies of the graph 
Dij and then identify nodes a^'s, bj q 's with the corresponding nodes in Ai and Bj (in G). 
Note that the A copies are sharing the same set of nodes, Ai and Bj, but other nodes are 
disjoint. 

3. Let G = (V, E) be the obtained graph (See Figure 5 for an illustration). 

The analysis: The next lemma shows that the size of an optimal solution in Ground PDS is exactly 
one more than the size of an optimal solution in the MinRep instance. The number of nodes in the 
constructed graph is at most I V(G r ) | < 4+ | V(C?)| + 10A \E(G)\. This shows that the above reduction 
is a gap preserving reduction from MinRep to Ground PDS with the same gap (hardness ratio) as the 
MinRep problem. Therefore the following lemma will complete the proof of the above theorem. As 
we mentioned above, the reduction given here is similar to the one used for proving the hardness of the 
directed PDS problem [13]. One important part of the above construction (see Figure 4) is the gadget 
on the set of nodes {a, (3, 7}. Note that there should be such a gadget between the center node in Dij 
and each u q and v q , in Figure 4 not all of the gadgets are shown (for example between u 2 and the 
center node). This gadget introduces direction into undirected construction, it allows the propagation 
in only one direction. After the center node of Dij is power dominated, all of the other nodes in Dij 
get power dominated (by the propagation rule). On the other hand, the power domination cannot 
propagate through the gadget in the other direction (toward the center node). 
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Figure 4: The graph D 



i.i 



Lemma 5.3 The pair {A*,B*) is an optimal solution to the instance G = (A,B,E) of the MinRep 
problem if and only if TI* = A*L)B*U{w*} C V(G) is an optimal solution to the instance G of £ -round 
PDS (for all I > A). 

Proof: The node w* should be in any optimal solution in order to power dominate , w;^ , , since 
otherwise we need to have at least 2 nodes from the set {wl, w^, w^} to get a feasible solution. By 
picking w* , all the nodes in A U B (and also the nodes inside -Djj's that are the neighbors of w*) will 
be power dominated. 

Assume that A* U B* is an optimal solution for the MinRep instance G. Now let us show that 
II = A* U B* U {w*} is a feasible solution to the PDS instance G. As described above, the nodes in 
AL)B and some nodes inside D^s are power dominated by w*. Consider a super edge AiBj in TL. Note 





(a) MinRep Instance G (b) ^-round PDS instance G 

Figure 5: The hardness construction 
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that the set A* U B* covers all the super edges in TC. So there exists an edge e q = {a-i q , bj q } E E(G) 
such that cn q E A* and bj E £?*, which covers the super edge AiBj. Since ai q and are in the set 
II they will power dominate their neighbors, u q and v q , in all of the 4 copies of D^ in G. After u q 
and v q are power dominated, the node d q will power dominate the center node in D^. It is easy to 
check that after the center node is power dominated, all of the nodes in Dij will be power dominated 
(through the gadgets on the nodes a,/3, 7). This shows that II is a feasible solution for PDS in G. 
Also it is straightforward to check that such a solution will power dominate the entire graph in at 
most 4 parallel rounds. Therefore, Opt £ (G) < \A* U B*\ + 1. 

By Proposition 2.3 the size of an optimal solution for PDS is a lower bound on the size of an optimal 
solution for Ground PDS. So it is enough to prove that the above upper bound is also a lower bound 
for PDS in G. Let II* C V(G) be an optimal solution for PDS. As we saw above, w* should be in any 
optimal solution for PDS. Now define A' = A n II* and B' = B n II*. First we prove that any optimal 
solution of PDS only contains nodes from A U B U {w*}, and then we show that (A 1 , B') covers all the 
super edges. Suppose for the contradiction that IT* contains some nodes which are not in AL)BU{w*}. 
So there are some D^s that cannot be power dominated completely by II* n (A U B U {w*}). By 
symmetry all the 4 copies of Dij are not completely power dominated. So the optimal solution II* 
needs to have at least one node from at least 3 of the 4 copies, and the remaining one might be power 
dominated by applying the propagation rule. By removing these 3 nodes from II* and adding ai q E A4 
and bj E Bj to IT* for some arbitrary edge e q = {cii , bj q } E E(G) we can power dominate all of the 4 
copies of . This is a contradiction with the optimality of IT* . This proves that any optimal solution 
will consist of nodes only from iUBU {w*}- To show that (A', B') covers all the super edges, it is 
enough to note the following: suppose no node from the inside of any copies of D^ is in the optimal 
solution; then any Dij can be power dominated only by taking both end points of an edge between 
the corresponding partitions (Ai,Bj). This shows that the size of an optimal solution for PDS on G 
is at least the size of an optimal solution for the MinRep problem on G plus 1. This completes the 
proof of the lemma. □ 

6 Integer programming formulations for ^-round PDS 

In this section we present an integer programming (IP) formulation for the ^-round PDS problem, and 
then we present a related integer programming formulation for the original PDS problem. Finally, we 
consider LP relaxations of these two IPs, and we show that they both have integrality gap of Q(n). 

6.1 An integer programming for £-round PDS 

Here we consider an IP formulation for the Ground PDS problem. Given an undirected graph G = 
(V,E) and a parameter 1 < £ < n, where n = \V\, define the set of parallel rounds T = {1, ...,£}. 
The variables in the IP formulation are as follows. Let S* be the optimal solution. We have a binary 
variable x v for each node v, that is equal to 1 if and only if the node v is in S* (S* = {v E V : x v = 1}). 
For each node v and a parallel round t £ T we have a binary variable z l v , where 2* = 1 means that 
the node v is power dominated on (or before) the parallel round t. For each edge {u, v} E E and a 
parallel round t E T we have binary variables Y*^ v and Y*^ u , where Y*^ v = 1 means that u can 
power dominate v at the parallel round t + 1. Before stating the IP formulation formally, we describe 
the constraints informally. We have a set of constraints for the termination condition saying that every 
node should be power dominated at the end of the last round (see (1)). The second set of constraints 
are for power dominating all the nodes in the closed neighborhood of a node that is in the optimal 
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solution (see (2)). We have another set of constraints for each edge {u,v} E E that checks if the 
propagation rule (R2) can be applied to u in order to power dominate v in the next round; node u is 
ready to power dominate v only if u and all of its neighbors except v are already power dominated 
(see (3)). The last set of constraints are for checking if node v is power dominated at time t; node v is 
power dominated only if either it is in the optimal solution or at least one of its neighbors can power 
dominate v at time t — 1 (see (4)). Once v is power dominated, it should remain power dominated. 
The term -\-x v in the right hand side of constraint (4) is needed to ensure that the z variables are 
monotone, that is, > z\. 





min x v 

V 




s.t. 






(1) 


i<4 


E V 


(2) 




Vu G V 




u£N[v} 




(3) 


Y l < z l 


V(«, u) : {«, v} e Vu; E JV[u] \ {u} , Vt E T 


(4) 


z l < V y*" 1 + x 


E V,Vt E T \ {1} 




uGN(v) 




(5) 


all variables are binary 





To get an LP relaxation, we relax the variables to be non-negative instead of being binary. It turns 
out that this LP relaxation is very weak. We can add the following valid constraints that make the 
LP stronger 2 . The first one forces the number of power dominated nodes at the first round to be at 
least equal to the size of the smallest closed neighborhood, and the second one forces the number of 
power dominated nodes to increase by at least one at each round. 

(6) 5(G) + 1<J24 

V 

(7) ^4 _1 + 1<E^ VtET\{l} 

V V 

We now prove that the new LP relaxation has integrality gap of £l(n). 

Theorem 6.1 Let £ be a given parameter, where I = O(logn). Then the LP relaxation for the l-round 
PDS problem has an integrality gap o/0(n) even on planar graphs. 

Proof: Consider the graph G that is obtained from the cycle on m nodes, C m , by creating a new node 
v' for each node v and connecting it by an edge to the original node (adding the edge {v, v'} to G) 
(see Figure 6). The graph G has n = 2m nodes and 2m edges. It follows from the proof of Theorem 
4.1 that the size of a minimum power dominating set in G is at least and the optimal solution 
power dominates the graph G in two parallel rounds. Therefore the size of any optimal solution for 

2 These constraints increase the optimum value of the LP on the cycle with 9 nodes, Cg, for £ = 3 from 0.6 to 1. 
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Ground PDS is at least |~^] for any I > 2. Now we show that the LP relaxation has the optimum 
value of O(l), and this completes the proof. Let U denote the set of nodes in C m in the graph G, 
and V denote the set of nodes of degree 1 in G (the newly introduced nodes). We assign value a > 
to all of the variables corresponding to the nodes in U, and value (3 > to all of the other variables; 
Vu € U : x u = a and £ V : x v = (3. Before giving a feasible solution with the objective value of 
O(l) (for the LP based on a, f3), we compute the value of the variables based on a and (3 for the first 
few rounds. 

If we apply the set of constraints (2) to all the nodes we get: Vu € U : z\ < 3a + (3 and Vt> G V : 
z\ < a + f3. Let u, v, w be the nodes of G as shown in Figure 6. By symmetry, we get 3 different 
types of Y variables: Y* r _ yu , Y*_> u and Y^^ u . Again the symmetry of the graph and the constraints 
imply that the value of Y only depends on its type. It is easy to check that the set of constraints (3) 
give the following inequalities: Y^ v < 3a + (3(= z*), Y„_> u < a + (3(= zl) and Y^ u < a + (3(= zl). 
Next we apply constraint (4) and we get: z\ < (a + (3) + 2 x (a + f3) + a = 4a + 3[3{= 3zl + a) and 
z% < 3a + j3 + /3(= z} l + f3). We can continue in this way, until we reach the parallel round t = I. 
It is easy to prove by induction that the following assignment for Y and z variables satisfy the set 
of constraints (2) to (4): Y*_^ v = z w Y^_ tu = = zl, where z u and z v are defined recursively as 

follows: 

zl +1 =a + 3zl (El) 
4 +1 =/3 + 4 (E2) 

Note that we assign the same value, 2*, to all nodes in U, and the same value, zl, to all nodes in 
V. Also we assign the same value for all Y variables of the same type. By combining the above two 
recursive equations El and E2, we get the following independent recursive definition for z u . 

4+ 2 = a + 3(3 + 3zl (E3) 
zl = 3a + (3 (E4) 
z 2 u = 4a + 3/3 (E5) 

It remains to check if the above assignments also satisfy the set of constraints (6) and (7). 

1 2 

(6) : m ■ (3a + (3) + m ■ (a + (3) > 2 => m ■ (4a + 2/3) > 2 2a + (3 > — = - 

m n 

1 2 

(7) : m ■ z l ~ l + m ■ zt~ l + 1 < m ■ z l ~ l + m ■ (3 + m ■ 3zt~ l + m- a=^a + /?> — = — 

m n 

Therefore, in order to satisfy the set of constraints (6) and (7), take a, (3 > be real numbers such 
that a + (3 = — . Equation E3 for any t > 1 implies that z^ +2 > 3 x z l u . By solving this recursive 
inequality we get z^ k > z^ ■ 3 k > ^ x 3 k . Hence, for k = [~log 3 m] we get z^ > 1, and consequently 
we get Zy k+l > z^ k > 1 by Equation E2. Therefor by taking I = 2 • |~log 3 m] + 1, the set of constraints 
(1) is also satisfied. Note that the above feasible solution has the objective value of m x (a + (3) = 1, 
since \U\ = \V\ = m. This shows that the LP has the optimum value of 0(1) and this completes the 
proof. □ 

6.2 An integer programming formulation for PDS 

As we mentioned before, the Ground PDS problem for £ > n — 1 is the PDS problem. Here we 
consider an integer programming (IP) formulation for the PDS problem that is different from the IP 
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Figure 6: A planar graph with large integrality gap 



formulation considered in the previous section for Ground PDS. The IP formulation here is based 
on finding an ordering in which the nodes are power dominated in the optimal solution. Given an 
undirected graph G = (V, E) with n nodes, let T = {1, 2, . . . , n}. In this IP, all of the variables have 
the same definition and meaning as before except z v . Here, z v indicates the round in which v is power 
dominated. The variable z v is equal to 1 if and only if v is power dominated at round t. Before stating 
the IP formulation for the PDS problem formally, we describe the set of constraints informally. There 
are two sets of constraints for presenting the ordering in which the nodes are power dominated (see 
(1) and (2)). There is another set of constraints for checking if node u can power dominate node v at 
round t. Node u can power dominate v only if either u is in the optimal solution or all of the nodes 
in its closed neighborhood except v are already power dominated (see (3)). The last set of constraints 
is to check if a node can be power dominated at round t. Node v can be power dominated at time t 
only if either it is in the optimal solution or at least one of its neighbors was ready to power dominate 
it at the end of the previous round. 

(IPo) min ^ x v 

V 

S.t. 

a) E^ = 1 

t 

(2) E^ = 1 

V 

t 

(3) Y u ^ v < z w + x u 

t'=l 

(4) 4< Yl Y uZl+x v 

u£N(v) 

(5) All variables are binary 



W € V[G) 

v< e t 

V(u,«) : {u,v} G E(G),Vt eT\{)i},Vw G N[u] - v 

\/veV(G),Vt gT\{1} 



It is easy to verify that the above IP formulates the PDS problem. Now we consider the LP relaxation 
that is obtained by relaxing the integrality of variables to nonnegativity constraints. We add the 
following valid inequality to the above IP to get a stronger LP. This new inequality forces a power 
domination step to occur at each round. 

(6) 1< £ YU V VtGT\W 
{u,v}eE 
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Now we show that this LP relaxation has a big integrality gap. 
Theorem 6.2 The LP relaxation for PDS has an integrality gap of Q(n). 

Proof: Consider the graph G = (V, E) that is obtained from the cycle on m nodes, C m , by attaching 
a node of degree one to each node of the cycle (see Figure 6). Clearly, G has n = 2 x m nodes. It is 
easy to check that the size of an optimal solution for PDS on G is at least y = § = Q(n). 

Now we show that the LP relaxation has an optimum value of 0(1). The graph G has n nodes 
and n edges. Assign value - to each variable x v and z l v and value ^ to each variable Y^ v . It is 
easy to check that all of the constraints are satisfied by this assignment. The objective value of this 
assignment is n x — = 1, and this implies that the optimum value is 0(1). □ 
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Appendix 



A.l Proofs 

Proof of Proposition 4.1: It is known that Dominating Set is an NP-hard problem on planar 
graphs [18], so we just need to prove the theorem for i > 2. We use (almost) the same reduction 
that has been used to prove the NP-hardness of PDS on planar graphs [9, 10]. Given a planar graph 
G = (V,E) and a parameter £ > 2 construct a graph G' , an instance of Ground PDS, as follows: for 
each node v E V create a path P v of length i — 1 and identify an end point of P v with the node v. 
This transformation preserves the planarity, and the theorem follows from the fact that the size of an 
optimal solution for the Dominating Set problem on G is equal to the size of an optimal solution 
for the Ground PDS problem on G' . The proof goes in the same way as the proof for the NP-hardness 
of PDS by Kenis et al. [9] and Guo et al. [10]. 

Assume that S is an optimal solution for the Dominating Set problem in G. Now it is easy to 
see that S power dominates G' in exactly I parallel rounds; the set S in one parallel round power 
dominates all nodes in V and then in the remaining I — 1 parallel rounds each node v E V starts to 
power dominate the nodes on its attached path, P v . 

Now assume that S' is an optimal solution for the Ground PDS problem on G ' . It is easy to see 
that there is such an S' that contains only nodes of degree at least 3, so consider that S' . Assume 
that S' is not a dominating set for G, so there is a node v E V that is not power dominated in the 
first parallel round. This means that v is power dominated by applying the propagation rule to one 
of its neighbors say u. It is not hard to see that u E V. Let v! be the neighbor of u in the path P u . 
Since u is not in S' , v! should also be power dominated through u by applying the propagation rule. 
This means that u power dominates both v and v! , but this is impossible. Therefore S' should be a 
dominating set for G. □ 

In the reduction for the NP-hardness of PDS in [9, 10] only an edge, that is, a path of length 1, 
is attached to each node. Here, we attach a path of length t — 1 to make an instance such that the 
optimal solution needs exactly t parallel rounds to power dominate the entire graph. 

A. 2 Dynamic programming 

Our dynamic programming is based on valid timed-orientations which is an extension of the new 
formulation for PDS introduced by Guo et al. [10]. It is similar to the dynamic programming for the 
Dominating Set problem given in [25] and the dynamic programming for PDS given in [10]. 

Fix a parameter I and consider the £-round PDS problem. Assume that the graph G = (V, E) and 
V' C V and a nice tree decomposition {{Xi C.V \ i E 1} ,T = (I, F)) of G with tree-width k are given 
as input. Let denote the subtree of T rooted at node i E I, and Y. L denote the set (UjW(T;) Xj)\Xi. 
Furthermore, let Gj be the subgraph induced on Yj U Xi, i.e. G{ = G[Yi U Xj\. Also denote by G^ 
the subgraph induced on JQ. Let n% and m; be the number of nodes and the number of edges in 
G\ respectively. The dynamic programming works on the bottom- up fashion. On each bag Xi, it 
considers all valid timed-orientations of the subgraph Gi and stores the number of origins together 
with the orientation on the edges of Xi as the states of the bag X{. In the other words, the valid 
timed-orientation is stored through the states of the bag. 

The state of a bag: The state s for a bag Xi defines the orientation of the edges inside G\, the 
time label assigned to the nodes of Xi, the number of directed edges from v E Xi to all nodes in Yj, 
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and also the maximum of the time-label assigned to the neighbors of v in the set Yi. In a bag state 
s we denote the state of an edge e = {u,v} G E{G' i ) by s(e), the time-label assigned to v G Xi by 
st(v), the number of incoming edges from Yi to v by S-(v), the number of outgoing edges from v to 
Yi by s+(v) and the maximum of the time-label assigned to the neighbors of v in Yi by s y (v). Let 
e = {u,v} be an edge in G\, then s(e) takes one of the following 3 values: "u — > v" , u v — > u" or 
"_L"; where the first two values shows the direction of the edge e in the valid timed-orientation and 
the third one indicates that e is left undirected. Consider a node v € Xi, st(v) takes a value from 

{0} U {1, 2, . . . ,£} U i 1, 2, . . . ,i\ U {+oo}. The node v with s(^) = a means that we still ask for a 

propagation rule to be applied to a neighbor of v and power dominate it, and s(v ) = a shows that the 
node v is already power dominated at the current stage of the algorithm. We use ||s(v)|| to denote the 
integer value of the label s(v ) ignoring the hat notation (e.g. ||2|| = 2). Also S-(v ) takes a value from 
{0, 1}, s + (v ) takes a value from {0, 1, 2} and s y (v) takes a value from {0, 1, ...,£} U {+oo}. Values of 
or 1 for s_(-y) and shows the exact number of incoming/outgoing edges to/from v , but 2 means 

that there are at least 2 edges. Let us denote by Si the set of all possible states for the bag Xi. It is 
straightforward to check that the number of bag states for Xi is |«Sj| = 3 mi x (2£ + 2) ni x 5™ 1 x (£ + 2) ni . 
Note that a node cannot have S-(v) = 1 and = 2 at the same time (this follows easily from 

Definition 3.4), so S-(v) and s+(v) have 5 different combinations. 

For each bag Xi we will compute and store a mapping Ai : Si — ► N U {+oo}. For a bag state s G Si, 
the value Ai(s) shows the minimum number of origins in the optimal valid timed-orientation of the 
subproblem induced on Gj under the restriction that the orientation of edges and labeling of nodes in 
Xi is defined by the state s. A bag state s G Si for the bag Xi is called invalid if 



(-■PI =) (3v G V n Xi : s t {v) = +oo) V 

(-.P2 =) (3v G Xi : (1 < \\s t (v)\\ < t) A (dr(u) + s_(«) > 1)) V 

(-P3 =) (3v G Xi : s t (u) = +oo A (dr(v) + s_(«) + dt(«) + s+ (v) > 1)) V 

(^P4 =) (3v € Xi : s t (v) = A d~(v) + 8 -(v) > l) V 

(-.P5 =) (3e = G E(GJ) : a(e) = "n -» u" A ((||at(«)|| = 1 A ||a t («)|| + 0) V (|| St (v)|| > 1 



where d~ (u), df(u), and A^[«] denote respectively the in-degree, out-degree, and closed neighborhood 
of u in the graph that is obtained from G[ by orienting edges according to the state s. Recall that 
the in-degree/out-degree shows the number of directed incoming/outgoing edges, but for the closed 
neighborhood we consider the undirected graph G^. Informally a bag state is invalid if it either violates 
any one of the valid timed-orientation's properties (PI to P5), or it cannot be extended to a valid 
timed-orientation. Now we describe our dynamic programming: 

Step 1: (Initialization) In this step for each leaf node i in the tree T, we define (initialize) the 
mapping Ai for each s G Si as follows: 



Step 2: (Bottom-Up Computation) In this step we compute in the bottom-up fashion from leaves 
to root the mapping corresponding to each bag in the tree. Recall that the tree nodes have three types: 



A||s t (u)|| < 1 + max {s y {u)} U {\\s(w)\\ : w G Ni[u] - v}))) V 
3v G Xi,3a G {i,... ,l\ : (s t (v) = a A d~(v) + s_(«) + 0)) V 
3a G {1, ...,£}'. s t {v) = a A d~ (v) + S-(v) = 0) 




\{v G Xi : s(v) = 0}| otherwise 



if either s is invalid or (3v G Xj : S-(v) + s+(v) + s y (v) ^ 0) 
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Join Node, Insert Node, Forget Node. In the following we describe how to compute Ai in each 
of these three cases. In each point of the algorithm we preserve the following invariant: for each s G Si 
there exists a valid timed-orientation for Gi which is compatible with s and has minimum number of 
origins equal to Ai(s) where the power domination of all nodes in Gi are justified except for the ones 
in {v G Xi : s(v) = a for some a G {1, . . . , £}}■ 

Forget Node: Suppose i is a forget node with child j, and assume that Xj = Xi U {x}. The bag 

p 

states s € Si and s G Sj are called forget- compatible and denoted by s ~ s , if 

(Fl) Ve G E(G'j) : a(e) = s'(e), 
(F2) W G : s t (v) = s' t (v) 

(F3) Vw G V(G'j) : s-(v) = s'_(v) + [s'{{x,v}) = "x -► «"] A s+(«) = s^O) + = "« -» x"] 

(F4) s'(cc) G {0,l,...,£}U{+oo} 



(F5) G y(G;.) : s y (t;) 



max{ S ;(t)), if {x, v} G S(G^) 

otherwise 



where [P] is equal to 1 if P is a true statement and otherwise. Now we compute the mapping Ai for 
the bag Xi as follows: Vs G Si 



Ms) 



+oo if s is invalid 

min{j4j(V) : s' G Sj, s ~ s'} otherwise 



Note that since x is not in Xi, by property (3) of the tree decomposition (Definition 3.1) it will never 
appears in any bag in the rest of the algorithm. This implies that the power domination of x should 
be justified within Gj. 

Insert Node: Suppose i is an insert node with child j, and assume that Xi = Xj U {x}. We introduce 
a mapping <f> : Si — > Sj in the following way. We map a given state s G Si to s' = cfi(s) G Sj in the 
following way. 



(11) Ve G E{G' j ) : s'(e) = s(e) 

(12) W G V{G'j) : s' t (v) : 



d if st(v) = a and [s({x,u}) = "x 
st(v) otherwise 



(13) W G V{G'j) : s'_{v) = s-(v),s' + (v) = s+(v),s' y (v) = s y (v) 
We now compute the mapping Ai for the bag Xi as follows: Vs G Si 



Ms) 



+oo if s is invalid or s_(x) + s+(x) + %(x) 7^ 

Aj((f>(s)) + [s(x) = 0] otherwise 



Again note that, since x appears in Xi but not in Xj by property (3) of tree decomposition it is the 
first time that it appears in the subtree rooted at node i. Also note that x cannot have a neighbor in 
Yi. 



Join Node: Suppose i is a join node with j and k as its children, and assume that Xi = Xj = Xj^. We 

say s' G Sj and 
conditions hold: 



say s' G Sj and s" G S^ are join-compatible with s G Sj and denote it by s ~ (s', s"), if the following 
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(Jl) Vv G V(G[) : (s.(v) = s'_{v) + s'L(v)) A (s+(v) = s' + {v) + A (s y {v) = max{sf y (v),a^(v)}) 

(J2) VveXi-. s t (v) g {0} u {i, . . . ,/J- U {+00} s t (v) = s' t (v) = s'/(v), 

(33) VveXr. s t (v) =a€{l,...,£}=> ((\\s' t (v)\\ = \\s'i(v)\\ = a) A (s' t (v) = a V s»(v) = a)) 

(34) Ve G E'(G-) : (s(e) = "±" s'(e) = s"(e) = "_L") A 
(s(e) = "u u" =► (s'(e) ^ "u -> u" A s"(e) ^ "u -» «")) 



Informally s is join-compatible with the pair of states (s\ s") if they are assigning the same label to 
the nodes and the same orientation to the edges, and also if a node is power dominated it should be 
justified in either s' or s". We now compute the mapping Ai for the bag Xi as follows: Vs G Si 



Ms) 



+00 s is invalid 

min ^Aj(s') + A k (s") - \{v G Xi : s(v) = 0}| : s' G <Sj, s" G s ~ (s', s")} otherwise 



Step 3: (At root r) Let r be the root of the tree decomposition T. Finally we compute the number 
of origins in the optimal solution for Ground PDS on the instance (G, V) in the following way: 

min {A r (s) : s G S r ,Wv G X r : s(v) G {0, 1, ...,£} U {+00}}. 



Recall that the number of bag states is \Si\ = 3 mi x (21 + 2) n ' x 5 n * x (£ + 2) n \ Let m e be the 

maximum number of edges that a bag can have, and also note that the maximum number of nodes 

is (k + 1). It is easy to check that each step of the dynamic programming can be computed in time 

0(c m<!+fclog ^), for some global constant c. This shows that the total running time of our algorithm is 
Q^ c m e +kloge . \y\\ 
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